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Materials exhibiting negative differential resistance have important applications in technologies involving mi¬ 
crowave generation, which range from motion sensing to radio astronomy. Despite their usefulness, there has 
been few physical mechanisms giving rise to materials with such properties, i.e. GaAs employed in the Gunn 
diode. In this work, we show that negative differential resistance also generically arise in Dirac ring systems, an 
example of which has been experimentally observed in the surface states of Topological Insulators. This novel 
realization of negative differential resistance is based on a completely different physical mechanism from that of 
the Gunn effect, relying on the characteristic non-monotonicity of the response curve that remains robust in the 
presence of nonzero temperature, chemical potential, mass gap and impurity scattering. As such, it opens up 
new possibilities for engineering applications, such as frequency upconversion devices which are highly sought 
for terahertz signal generation. Our results may be tested with thin films of Bi 2 Se 3 Topological Insulators, and 
are expected to hold qualitatively even in the absence of a strictly linear Dirac dispersion, as will be the case in 
more generic samples of Bi 2 Se 3 and other materials with topologically nontrivial Fermi sea regions. 
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1 Introduction 

Topological Insulators (TIs) are a new class of materials with a fully insulating gap in the bulk but gapless (conducting) 
Dirac fermion states on the surface 1-3 , and have garnered tremendous interest in condensed-matter physics, material 
science and electrical engineering communities 1-14 . Recent experimental realizations of TI states in compounds like 
HgTe, Bii-^Sb^ and Bi 2 Te 3 fueled further enthusiasm in their possible application in devices. The Dirac cones 
on the surface of 3-dimensional TIs are reminiscent of the Dirac cones in 2-dimensional Graphene, another exotic 
material which has attracted considerable attention 15-21 . Notably, Graphene has been theoretically predicted 19-20 
and subsequently experimentally shown 21 to exhibit strong nonlinear electromagnetic response owing to its unique 
linear Dirac dispersion 15-21 . 

Inspired by these properties of Graphene, we ask if similar, if not more desirable, nonlinear behavior is also 
present in the TIs. As we will show in this work, the answer is in the affirmative: In fact, the response curve of TIs 
is even more nonlinear, with an exotic regime of negative differential resistance that persists even in the absence of a 
strictly linear dispersion. More precisely, the response curve takes the form of an ’N’ shape with negative differential 
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resistance in the middle segment, similar to the shape of the response curve leading to the Gunn effect in GaAs 22 , 
but due to a completely different physical mechanism. This suggests an array of potential optoelectronics applications 
beyond those of Graphene and the Gunn diode. 

The enhanced nonlinearity of the response of TI surface states can be understood as follows. A TI heterostruc¬ 
ture has two conducting surfaces, the top and bottom surfaces, while a Graphene sheet only has one. Due to the 
existence of the substrate in the TI heterostructure, structural inversion symmetry has to be broken, leading to a break¬ 
ing of the degeneracy of the two TI surface states which opens up a Rashba-type spin splitting 23 . This results in an 
unique Dirac ring bandstructure which exhibits a much stronger nonlinear electromagnetic response than a Dirac cone 
alone, thereby opening up a venue for interesting physics as well as potential applications. 

In this work, we model a TI heterostructure as a Dirac ring system, and analytically and numerically study its 
semiclassical nonlinear response. We first consider the case of an ideal Dirac ring, i.e. at zero mass gap, temperature 
and impurities. Next we consider deviations from these ideal conditions, and crucially show that the characteristic 
features of the response curve remain robust. We further discuss how this semiclassical analysis can be generalized to 
a more general setting with scattering and/or Rashba-like dispersion, and its implications for the output signal. Finally, 
we discuss some experimental proposals and engineering applications. 


2 Results 


We first introduce some basic theory on the semi-classical electromagnetic response of a generic Hamiltonian. With 
that, we present our main results on the characteristic nonlinear response curve of Dirac ring systems. Such systems 
have been detected in the surface states of thin films of B^Ses TIs via ARPES experiments 23 , and we will return to 
discussing the experimental signatures of our results after developing its theory. 


Theory of semiclassical response Consider a generic system described by a Hamiltonian H (p) under the influence 
of a driving field E(t). At the semi-classical level, the field shifts the crystal momenta p of the partially occupied 
bands, leading to an induced current 


J(t) = e{v) t = 


( 27 rh)‘ 


j d 2 pg{p, t)i 


( 1 ) 


where v = V^e(p) is the canonical velocity and e(p) is the eigenenergy for a particular band. In other words, J is 
the expectation value of the current ev over states weighted by the time-dependent equilibrium occupation function 
g(p,t). As will be shown rigorously in Sect. |4j the equilibrium occupation function takes the form of the Fermi- 
Dirac distribution g(p,t) = F 0 (p — Po(t )) 9 but with momentum shifted by an effective driving impulse po(t) = 
e f l E e ff(t')dt'. We shall elaborate on exactly how E e ff depends on E and the scattering time r in Sect. 


2.0.5 


For 


now, we shall proceed by treating p 0 as an external influence, and note that E e ff — E in the limit of zero scattering. 
This approach considers only intra-band transport processes, and is valid when po originates from an oscillatory 
electric field with h£l <C max{/i, T, m}, i.e. with frequencies under 20 THz in typical applications. Henceforth, we 
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shall work in units where 2nh = 1 and kp = 1 for notational simplicity. 


If the bands of the hamiltonian are isotropic with eigenenergies e(p) = e(p) where p = \p\, Fq depends only 
on p via e(p) and we can further express J as 

J = e J d 2 pF 0 {e(p-p 0 ))Vpe(p) 

= e J d 2 pF 0 (e(p))Vpe(p + p 0 ) 

= e J j(p) dp (2) 

where 

3 (?) =pFo(e(p)) J Vpe(p + p 0 )d9 (3) 

is the contribution from the p-momentum shell. This decomposition is particularly useful since j (p) can often be 
expressed in closed-form, although J usually cannot be. 

We now specialize to the Hamiltonians with a Dirac ring, which is the focus of this work. The energy dispersion 
of a Dirac ring system takes the form (Fig. [T} 

e(p) = ±^/m 2 + (vpp — M) 2 (4) 

with vp the Fermi velocity, m the (half) gap and ^ the ring radius. It is known as a Dirac ring system because it is 
rotationally invariant about 6 = tan -1 and has a ring of band minima with linear dispersion at radius p = 
Note that it reduces to the Dirac cone in Graphene when M = 0, which was systematically studied in Ref. [19]. In TI 
heterostructure realization reported and analyzed in Refs. [23-26], m is the mass gap induced by interlayer coupling 
and M measures the extent of inversion symmetry breaking. 


2.0.1 Nonlinear response of Ideal Dirac rings 


We shall first study the proptotypical case of ideal Dirac rings, where the gap m and temperature T are both zero. In 
the limit of small chemical potential p <C M, which can be tuned by varying the gate voltage 26 , the filled states on 
the valence band form a very thin ring bounded by inner and outer Fermi momenta pi = M ~J l and Pf = M +J L , he. 
with radial thickness 2 pb/vp- Being such a thin ring, its total current is thus well-approximated by j(M/vp) in Eq. [ 3 J 
which is simple enough to visualize schematically and exactly evaluate analytically (Eq. tut 

The canonical velocity is simply given by 

v(p)\m= 0 = Vp^O) = v F sgn(v F p - M)p (5) 

which is a vector of constant magnitude vp- It is always pointing in the radial direction, and is positive(negative) 
outside(inside) the Dirac ring. Such a scenario occurs when the filled states forms a non simply-connected ’ring-like’ 
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Figure 1: (Color online) Left) The bandstructure of a Dirac ring system ±e(p) given by Eq. [4j with asymptotic 
canonical velocity (slope) given by vp k- Right) 3D view of the same band structure, with the axis of rotation the energy 
axis and the plane of rotation the p x , p y plane. The occupied (red) band is separated from the partially occupied (blue) 
band by a mass gap 2m. They can be obtained by gapping out two intersecting Dirac cone with energy difference of 
2 M. A ring of filled states (red) at chemical potential p > m lies in the ring of minima on the valence (blue) band. 

shape instead of a ’blob-like’ shape. In this case, the lower dimensionality of the ring may allow for some novel kind 
of ’destructive interference’ to occur between the contributions on both sides of the ring, and hence lead to negative 
differential resistance. Before proceeding with the calculations, let us attempt to understand that more intuitively. 

Without a driving field, po = 0 and the ring of filled states lie exactly above the ring of Dirac nodes. On either 
side of it, the vector field v(p) points in equal and opposite directions, thereby resulting in a zero net current in Eq. 


reciprocal (momentum) space, leading to an imbalance between the contributions of (p) from inside and outside the 
ring. As illustrated in Fig. [2] J is relatively large for small po because its arises from v contributions that point in 
the same horizontal direction. As \po | increases to more than half the radius of the ring, contributions from inside the 
ring oppose those from outside the ring, thereby leading to a decrease in | J|. Finally, for p 0 larger than the radius of 
the ring, the filled states disentangle from the Dirac ring completely, and the v contributions point in same horizontal 
direction again, adding up more strongly than before. As po continue to increase, v will eventually become parallel, 
leading to a maximal current \J\ proportional to the size of the ring. 


33 andp] Upon a small impulse po = \po \ from the driving field, the ring of filled states will be slightly displaced in 
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Figure 2: (Color online) Left) For small driving impulse po , the current J = e(v) is relatively large as it is contributed 
by v whose horizontal components add constructively (light blue). Middle) As po increases to the range M < vpPo < 
2M, the horizontal components of v experience a partial cancellation (light blue), thereby lowering the resultant 
J = I J|. Right) Beyond po > 2 M/vp, the ring of filled states disentangles with the ring of Dirac nodes (grey). The 
contributions from v (light blue) add even more strongly, resulting in an even greater J. Very importantly, the behavior 
described here still holds for a ring thickened by (not too large) chemical potential or temperature, and even remain 
qualitatively true in the presence nonuniformities in \v\, i.e. from nonzero mass. 


The above arguments suggest a response curve J = \ J\ that rises sharply to a moderately large value when p 0 
is very small, decreases when M < vpPo < 2M, and increases to an even larger value for larger pq. We identify 
the decreasing region as the region of negative differential resistance °c < 0. This agrees exactly with analytic 
expression derived in section [4] and plotted in Fig. [3] 


Xn=o )Al <M = 4Q[(1 + Q) ( EllipticE [a] + EllipticE 


+(1 - Q) (EllipticF 


7T — (f 


— EllipticE 


7T + Lp 


7T — (p 


— EllipticF 


7T + P 


EllipticF 


l2' a l 


}/jep 0 


( 6 ) 


with Q = a = (l+Q) 2 ’ v = cos 1 2 ^, EllipticE^, A] = y/l - A sin 2 6d0 and EllipticF^, A] = -j=k==dd. 

Note that only the dimensionless combination Q = enters the expression. 
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Figure 3: (Color online) Plots of the current in Eq. [6] against vfPo/M (defined as Q~ x in Eq. 16), where po is 
magnitude of the effective driving impulse, po = eE 0 r in the DC scattering limit, where r is the scattering time. Left) 
Behavior of J/(evp) for p = 0.05, 0.1 and 0.2 in units of M. We see a strong linear regime at very small po, and a 
region of negative differential resistance ^ < 0 at larger values of po due to destructively interfering v. Right) The 
same plot normalized by p: J/(pevp) for p = 0.03, 0.1, 0.3 and 0.7 in units of M. We now clearly see the asymptotic 
properties of the varies curves. As p decreases, the response curve becomes ’sharper’ due to the decreasingly thickness 


of the ring of filled states, and eventually converges nonuniformly to G(Q, Q) (from Eq. 16) as p 0. 


2.0.2 Distortion of a sinusoidal signal 


Below, we show how an ideal Dirac ring system distorts periodic signals of different amplitudes. Due to the segment 
of negative differential resistance in the response curves in Fig. [3j additional kinks and lobes are introduced in the 
output signal J. These lead to larger high-frequency components than what can be obtained with Graphene 19 , whose 
distorted output signal is a square-wave which also appears in the A 0 oo limit of the Dirac ring (Fig. [4]). 
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Figure 4: (Color online) Output currents corresponding to input signals E e ff = Aq^^ cos fit, where A 0 = 
0.5,1 5 1*5,2,3,10 for very small fi. Note the almost square-wave-like waveform for small and large Ao, where the 
response curve is mostly flat (Fig. [3]). The negative differential resistance part of the response curve produces the 
interesting lobes in the output current. They will be smoothed out at larger fi, nonzero gap or nonzero temperature. In 
the Aq —>> oo limit, the Dirac ring becomes irrelevant and the response curve reduce to that of Graphene (oc sgn p 0 ) in 
the same limit. 


The extent of the lobes in the output can be quantified by the Fourier coefficients of the output current 

1 r 2n 

fn = ~ Jif) s in nt (7) 

ft Jo 

Assuming an input impulse oc sint, we find that only odd coefficients f^n-x are nonzero. The distribution of 
characterizes the frequency multiplication efficacy, which is an important concern in the production of low-power 
electromagnetic radiation from lower frequencies. This is of particular exigence in the frequency window of 0.3 to 20 
THz (commonly known as the Terahertz gap) where, despite a multitude of applications across engineering, material 
science and medical disciplines 19,27-29 , inexpensive and compact sources for the THz radiation are still lacking. 

It is interesting to compare the decay of for the Dirac ring with the slowest decay spectrum possible 

without negative differential resistance ^ < 0. The latter is given by a response curve proportional to sgnpo, 
which can be realized 19 in Graphene at low (or zero) chemical potential and temperature relative to vfPo- Simple 
computation yields a harmonic decay profile 2n - 1 > n = 1, 2,.... This is compared with that of the 

Dirac ring for various values of mass gap, temperature and input signal amplitudes in Fig. [5] We see a larger frequency 
multiplication factor, i.e. larger from the gapless Dirac ring across most of the harmonics / 3 , /s,..., / 19 . When 

an energy scale is introduced by either nonzero gap or temperature, / 2 n -i decays exponentially. Even then, the 
frequency multiplication factor still outperforms that of Graphene for the the first few harmonics. 
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Figure 5: (Color online) Ratios of the frequency multiplication factor of the Dirac ring with that of zero- 

temperature Graphene ( oc 2n 1 _ 1 ). The amplitudes A$ are given by E e ff = Aq^A cos Dt, which give maximal 
effective impulses p™ ax = A&M .. For reference, the continuous black line corresponds to Graphene at zero temperature 
T, which has the maximum possible frequency multiplication factor without negative differential resistance. It is 
exceeded by that of the gapless, zero-temperature Dirac ring at several different harmonics. This is expected from the 
appearance of additional, higher-frequency lobes in Fig. |4j where the sharper lobes in the Ao = 2 curve now manifest 
themselves as a suppression of the first higher harmonic relative to the A 0 = 1 case. At nonzero T or m, the 
higher harmonics are exponentially suppressed. However, the n = 2 orn = 3 multiplication factor can still exceed the 
maximal possible without negative differential resistance (black line), as well as that of Graphene at T = 300Ff, fi = 0 
under the same input signal as given in the discussion section. 


Response of imperfect Dirac ring materials So far, we have considered Dirac ring systems with vanishing mass 
gap m, such as those realized in B^Ses TI thin films heterostructures thicker than 6 quintuple layers (QLs) 23 . To 
demonstrate the robustness of the regime of negative differential resistance, we shall now analyze Dirac ring systems 
which possess a gap and correspondingly a departure from perfect linear dispersion. They arise in sufficiently thin TI 
films where hybridization of the surface states on either side of the TI opens up a gap due to wavefunction overlap 24-25 . 
This had been predicted 24-25 to occur and was indeed observed 23 in B^Ses heterostructures thinner than 6 QLs. 
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Real materials under experimental conditions furthermore experience non-negligible temperature effects and 
impurity scattering, both of which can undermine the preceding Dirac ring interference analysis. But very impor¬ 
tantly, we shall show that the qualitative features of the response curve, particularly the region of negative differential 
resistance, remain robust. As long as the occupied states still occupy a ring-like region in reciprocal space, we indeed 
observe: 


1. A rapidly increasing response J for small po ; 

2. A region of decreasing J < 0^ at moderate po « M/vp, the radius of the ring; 

3. An increasing, even larger J for larger po- 


To understand why, it is useful to examine Fig. [2] again. The non-monotonicity of the response (Center) is a con¬ 
sequence of the destructive interference of the contributions of v from inside and outside the ring. This is a generic 
feature for a ring of energy minima, since the v = Vpe(p), which is always of opposite sides of the ring. In par¬ 
ticular, note that it does not depend on the dispersion being linear, although the destructive interference will be less 
pronounced when the ring is thickened by large p or m, or fuzzied by nonzero temperature T. 

Below, we shall substantiate the above schematic arguments with detailed analyses and numerical results. 


2.0.3 Effect of nonzero band gap m 


In a gapped Dirac ring system, Eq. [4] gives the canonical velocity in a gapped Dirac ring system as 

vpp — M 


v(p) = Vpe(p) = 


\/m 2 + (vpp — M) 2 


vfP 


( 8 ) 


which vanishes linearly near the Dirac ring. As such, we expect a more ’rounded’ response curve, as shown in 
Fig. [6] (Left). Note that Eq. [8] also describes the small po response in a Rashba system. The latter, however, has 
a asymptotically quadratic dispersion which leads to a (rather uninteresting) linear response for large po > m/vp , 
which also suppresses the region of negative differential resistance. 


The results for nonzero mass gap m are shown in Fig. |6](Left). The response curves exhibit the same qualitative 
shape, but with a broader linear response regime for small , which is studied in more detail in Sect. [ 5 ] Examples 
of such cases include the bulk states of HgTe 31 and GaAs 32 quantum wells with inversion symmetry breaking. 
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Figure 6: (Color online) Left) The current response for different values of masses m, all at chemical potential /i = 
0.3 M and temperature T = 0. As expected, we observe the same qualitative, though ’rounder’, response curve with 
increasing m. The current diminishes to zero as m /i, when both \v\ and the thickness of the filled ring diminishes. 
Right) The current response for different temperatures, with M set to O.leV and /i = 0.3 M. There is essentially no 
difference in the curves at T = IK and T = 77 K (liquid nitrogen temperature). 


2.0.4 Effects of nonzero temperature T 


In the presence of nonzero temperature, the Dirac ring is smeared out over a radius of A p « T/vp. Furthermore, hole 
carrier from the —e(p) band also participate in current transport. The smearing of the Dirac ring results in qualitatively 
similar modifications to the response curve as increasing fi or m. However, the participation of the hole carriers 
undermines the destructive interference of v described in Fig. [2J and can destroy the characteristic non-mono tonicity 
and thus negative differential resistance of the response curve at sufficiently high T. Nonetheless, we numerically find 
(Fig. |6](Right)) that the ideal response curve retains its qualitative shape till room temperature (300K) for M ~ O.leV, 
the inversion symmetry breaking scale in experimental realizations. This cutoff temperature is approximately 3000M 
in units of K/eV . 


2.0.5 Effect of scattering 


Scattering is inevitable in real materials. For instance, a scattering time of r = 4.2 x 10 14 s was reported for a Bi- 
based TI compound 33 . Here, we shall model the effect of scattering by a classical relaxation time r, and show that its 
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effects can be incorporated into our calculations by replacing the input driving field E with a ’renormalized’ effective 
driving field E e f /. 


The electronic states are distributed according to a time-dependent occupation function g(p(t),t) that obeys the 
Boltzmann equation 


dg = dg dg .dp = fj^g 

dt dt dp dt t 


where / = f(p(t )) = F 0 (p(t)) = (1 + e^ >*)) 1 is the local equilibrium electron distribution, Fo being 
the Fermi-Dirac distribution which can deviate from g = g(p(t),t),\he true, non-equilibrium distribution. Being the 
equilibrium distribution, the functional form of f(p) remains unchanged, and all fluctuations in time occurs only 
through the time-varying momentum of an electron p(t). By contrast, g can deviate from this equilibrium through 
its explicit time dependence, although it will relax towards / with a characteristic time r. We will assume spatial 
homogeneity throughout. 


Eq. [9] has an explicit solution though elementary calculus: 

dg(t) _ C , , e (t '~ t)/r df(t') 
dt J t dt' 


( 10 ) 


where the dependence on p(t) has been made implicit. Intuitively, ^ is a Laplace transform of with contributions 
from earlier times exponentially suppressed due to scattering. As derived in Sect. [4] Eq. 10 has an explicit solution: 


g(p,t) = Fq i^p + e J E eff {t')dt^j ( 11 ) 

where 

Eef fit) = fdt'e^-^Ep. ( 12 ) 

is the effective electric field that contains the effect of scattering. Essentially, E e f / is contributed by past increments of 
the electric field E that are exponentially suppressed with a characteristic time r. As r increases, the increments have 
more time to constructively interfere. Physically, \E e ff\ increases as r increases because a longer ballistic motion 
contributes greater to the momentum shift of each particle. Of course, \E e ff(t)\ < ma x t ' <t \E(t')\. The manifest 
causality structure results from the time-reversal asymmetry of the system cf. Eq. [9] 


For a simple sinusoidal driving field E(t) = Eq sin fit, Eq. [l2](or Eq. [ 34 ]) can be solved for the exact effective 


field 


fir 


E e ff(t ) — Eq ^ ^ (cos fit H - Tfl sin fit') 


(13) 
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leading to the effective impulse 


P 


P~ 


= P~ 


erEp 
1 + (Or) 2 

erEp 


V^l + (Or) 


(sin Of — rO cos Of) 
sin(Of — tan -1 Or) 


In other words, scattering leads to a renormalization of the amplitude by a factor 


Qr 


(14) 


and a phase retardation 


corresponding to the duration 


tan Qt 


r [ 1 + 


l + (f2r 

^. The latter is trivial for a single frequency mode, but will lead 


n ' ^ 1 3 

to interference effects when there is a mixture of modes. That will further studied in Sect. |4j in the limits of weak and 


strong scattering. Note that an expression very similar to Eq. 14 also appeared in Ref. [19] in the context of radiation 
damping. But here r has a more generic interpretation, encompassing all generic damping mechanisms. 


The DC limit can be recovered by taking Clt — | <C 1, so that E(t) « Ep. The effective impulse is equal to 
eEpt in the ballistic regime, but saturates at eEpr when f is increased to r and effect of scattering is felt. 


3 Discussion 

Through an approach based on the Boltzmann equation, we have analytically and numerically obtained the novel non¬ 
linear response curve of a Dirac ring system representing Topological Insulator thin films. Due to the special topology 
of a ring-shaped region of filled states, the response curve is intrinisically non-monotonic, leading to a regime with 
negative differential resistance which heightens frequency upconversion of a periodic signal. This property remains 
robust in real, experimentally fabricated samples in the terahertz gap, with nonzero mass gap, chemical potential, 
temperature and impurity scattering rate. 

One test of our results will be to reproduce the output signals in Fig. [4] To produce the Ap = 1 curve, for 
instance, we can input either an electrical or optical sinusoidal signal 34 of amplitude Ep = ^^Ap on a 

sample of Bi 2 Se 3 thin film thicker than 6 QLs. Experimentally obtained parameters M « O.leV, vp = 4.52 x 
10 5 ra/s, 23 and r = 4.2 x 10 -14 s, 33 correspond to a realistic electric field of magnitude 34 Ep = 5.27 x 10 4 V/cm 
and frequency 35 Q = 100 GHz. The output current should closely agree with that in Fig. [5] at the temperature of 
liquid nitrogen T = 77 K, and deviate slightly from it at room temperature T = 300if (cf. Fig. |6|Right)). All 
other output signal shapes in Fig. [4] may be reproduced by varying the input signal amplitude. Since the effective 
impulse from the field increases linearly with the scattering time r, TIs with much larger r will exhibit the aforesaid 
negative differential resistance and nonlinear response behavior at much weaker physical electric fields. This will be 
very attractive for low power applications, and is well on track to becoming a reality with the rapid development of TI 
material growth techniques as well as the discovery of novel TI materials. 


4 Methods 
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The current response of a gapless Dirac ring Here we derive the analytic expression for the response of an (ideal) 
gapless Dirac ring, for which the p —» 0 result in Eq. [6]is a special case. 


Substituting the canonical velocity v(p) | m=0 = V^e(p) = v F sgn(v F p — M)p in Eq. [3j where Fo(e) = 0(—e) 
at temperature T = 0, we obtain 


j(p)\m=o = v F pF 0 (e(p)) /sg- {v F \p + p 0 \ - M ) 


P + Po 

Ip+P ol 


dd 


Pov F F 0 (e(p)) /Sgn( + £>o + 2/ppo cos 0 — M /'o 


p(p 0 + P COS 0) 
\/Po +P 2 + 2 PPo cos6» 


dO 


= VF F 0 (e(p))G(^,—) 
\Po v F po) 


Po 


(15) 


where G(P, Q), P = Q = is a closed-form expression given by 


G(P,Q) 

2 P 


= (1 + P) EllipticE [a] + EllipticE 


7T — (f 




+(1 - P) EllipticF 


7T — (/9 


— EllipticF 


— EllipticE 

7T + (/? 


7T + <p 




EllipticF 


L2 ,a 


(16) 


with a = p = Re 


i+p 2 -q 2 

2P 


, EllipticE[0, A] = y/l — X sin 2 #<7# and EllipticF[0, A] = , 1 =d0. 

u U V 1 — A sin 2 0 

The functional dependence G(P, Q) = G implies that of the 3 parameters p,po and M/v F , only two can 

affect the result independently. For instance, the effect of letting M 2 M is equivalent to that of po —>> po/2, 


p -F p/2 and G —>► 2G. 


We are now ready to calculate J proper. It is just a sum over annuli of different radii p , each contributing 
a current j (p). From Eq. |4j p defines two Fermi surfaces with inner and outer Fermi momenta pi = M ~J l and 
Pf = M +J x • Considering the interesting case of small p < the occupied momenta lie in a ring of inner and outer 
radii pi and p F . In the limit of p 0, the occupied states form a very thin annulus and an illuminating closed-form 
expression exists for J : 


J m= 


=o,/x<m — e 


L 


M + p. 
v F 


M-n 

F 


j(,p)\m=odp 2iiep 0 v F G 


M M \ 


; 


Po = 2peG(Q,Q)po (17) 


K v F p o v F p 0/ 

which is just Eq. [6] For larger p, numerical integration yields the plots in Fig. [3] which has the following properties: 


• The current is notably proportional to v F . This property is unique to the quasi-ID shape of the ring. While the 
velocity operator v is proportional to v F , the area of the ring is independent of it. This is because its radius is 
proportional to it, while its thickness is inversely proportional to it. 

l lf p > M, the occupied momenta lie in a disk of radius p F for the lower band, and of radius \pj | for the upper band. Note that all incompletely 

filled bands must be included. Of course, the model must be lattice regularized if we want to study huge p F of the order of i r. 
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• The response at any finite // does not converge uniformly to Eq. [6j At small perturbations po, the ring feels 
opposite velocity fields at both sides of the Fermi surface (FS), and there is a resultant linear regime. From Eq. 
[6j one can show that 

Jm=o^<M,p 0 <^ ~ 2ireMpo (18) 

for 0 < po < m t0 a ver y high degree of accuracy, as shown in Fig. [3] Next comes a nonlinear regime where 
the response is negative approximately in the range fi < 0.6^. 


At po = — or Q = 1, we pass a special point where the ring untangles from the ring of Dirac nodes. Here, 




m=0,ii<M,po 


27repM 

v F 


(19) 


• For larger p 0 or Q 1, we have 


It can also be shown that 


J, 


m=0,/z<M 


J m=0,/x>M — 


AnepM 
v F 

27re(/x 2 + M 2 ) 
VF 


( 20 ) 

( 21 ) 


Current response of a gapped Dirac ring Here we derive some analytic results for the response due to small per¬ 
turbations about a thin gapped Dirac ring, so as to understand how the ring structure affect the linear response of an 
otherwise massive system. 

For sufficiently small p — m and po, the inequality m \v F \p + Po | — M\ holds and the canonical velocity is 
approximately 

v(p) ^ —(v F p — Mp) (22) 

m 

This is also what we have in the vicinity of the ring of minima of a system with Rashba splitting. Substituting it into 
Eq. [3]like before, we obtain 

3 ip)\m large « (2t jpv F - MG(p/p 0 , 0)) p 0 (23) 

m 


The first term is linear in po, agreeing with that of usual materials with quadratic dispersion cx yp 2 ; upon doing 


the p integral, we recover 


r _ e7r v fP 2 f * 

J quadratic PO 

m 


(24) 


where p F is the (outer) Fermi momentum. The second term is a Dirac-like contribution that modifies the response 
from the lowest-order quadratic approximation. Integrating both terms, we obtain 


where p F , Pi = 
Eq.|3|and[T6| 


J | m large 
M±yjp?-rn 1 


evp 

m 


-(P F ~ P 2 l)P0 ~ — ( JDirac{PF ) JDirac(Pl)^j 


(25) 


and JDiracip) = ^PqV F f Q G(P , 0 )dP po denotes the Dirac cone (M = 0) current in 
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The effective driving field due to scattering Here we show that the non-equilibrium state distribution g is of the 
form 

g(p,t) = Fo^p + e j E eff (t')dt'^j , (26) 

and derive E e ff in terms of the original driving field E. The Boltzmann Eq. [9] has an explicit solution Eq. |To 1 


dg{t) 

dt 


= t dt' 

7 —OO 


o{t'-t)/r 


df(t ') 


dt' 


(27) 


where / = f(p(t)) = Fo(p(t)). To motivate the solution to Eq. 10 we first solve it for a periodic driving field. For 
each fourier component £2, we have = — eE{t ) = —eE^e 1 ^, so that 


dg(t) 

dt 


f 


rt df de dp 

/ —OO 

/ ,,/ c ' " u ’±Oy C yy J/ 

-e / dt ---u(p) • B^e 

J-oo t de 

dF 0 (e(p )) _. - /•* 

-e \^ J, v(p)-En dt'- 
7—00 


r de dp dt 

dFo(e(j>)) 


iQt' 


,(t' — t)/r+iQt' 


= —e 




de 

dFojejp)) 
de 

dF 0 (e(p)) 


v(p) ■ E n e 


vip) • Bn: 


/ 

7o 


00 -(i«+l/ r )T 

dT-- 

0 T 


a iVtt 


= —e 


de 1 + 

dFo(e(p)) v(p) ■ E(t ) 


de 


1 + i£2r 


(28) 


For a generic periodic driving field, Eq. 28 still holds, but with a sum over all fourier modes £2. Hence the effect of the 
r damping is the replacement of the fourier coefficients (Ed)n = which can also be guessed from elementary 

considerations. 


As such, let’s define a damped momentum pd which responds to the damped electric field Ed with fourier 
coefficients ( Ed)n • From Eq. 28 we obviously hav^] 


dg(p(i),t) _ df(p d (t)) 


dt 


dt 


-eV/ • E d 


This integrates to 


g(p, t ) = B 0 (V+ e J* (Ed(t') - E(t'))dt^j = F 0 (p + e J ‘ E eff (t')dt' 


(29) 


(30) 


2 Care has to be taken in handling these PDEs: while g depends on t both explicitly and implicitly through p(t), the dependence of / on t is only 
implicit through p and pd, respectively before and after integrating out the damping effect. 
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where 




o(t'-t)/r 


— / dt'Sfi') / -i— 

27r -/ J 1 + sb 

f m 

J* Me# 


-dt' - E{t) 


dEtf) 

dt 


4.0.6 The diffusion limit of small r 


( 31 ) 


In the limit of strong scattering, only recent memories of A E survive, and an expansion about t' = t in Eq. 


31 


dE(t) 2 d?E(t) o d 3 E(t) 
—t — - r ——-r 


dt 


— -r^ E(t) 

1 


d 2 £ 


(ft 3 


gives 


(32) 


The above expansion is valid in the regime r < At, where A t ~ ft, 1 is the characteristic time scale at which E 
varies. To linear order, the current is thus 


J linear ^ d pv(p)Fft C E e ff(t )dt 


e Jd 2 p ( F 0 (p) + e 


dMp) , 

de^ 


9 / j2 dp o (P) 


I 


d p- 


dtf, 


v ■ tE( t) + 


f E eff (t')dt' 
2 dE(t) ^ 


dt 


+ - ]v 


(33) 


This linear approximation coincides exactly with the usual derivation of the Drude formula. If we allow for large 
perturbations, we will instead have 


J non—linear 


i J d 2 pv(p)F 0 ( p + erE(t) 


+ er‘ 


,dE(t) 

dt 


(34) 


Clearly, the argument of Fo is the unperturbed crystal momentum p plus the impulse from E over an effective 
duration of r. Note that it is not a Taylor expansion of E(t — r); it is the series form of the exponentially suppressed 
field given in Eq. 31 with terms given by r n d dtn -i^ • 
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4.0.7 The ballistic limit of large r 


First, we consider what happens when r —>■ oo, which is the limit studied in Ref. [19]. Since the system is driven 
periodically, E(t') has equal positive and negative contributions, and 


lim 

T—>- OO 


dg{t) 

dt 



(35) 


That f = 0 forces the g to have the functional form F 0 (p — eA(t )), where explicit time dependence only enters 
through A(t) = — J t E(t')dt'. Fq is fixed to be the Fermi-Dirac distribution by considering the limit \E\ —>> 0. In a 
sense, this argument is a simple yet insightful semiclassical justification of the minimal substitution p p — eA for 
a distribution: only through this substitution will g remain a constant as we follow a particle, as it should be in the 
absence of any other force. 


Now, let us consider first-order scattering contributions ~ ^. From the second line of Eq. 

g(p,t) « F 0 (V- e J E{t')d£ + ^J* £(i')e (t '-‘ )/T ) 


31 


we have 


(36) 


The ^ term keeps track of the effects of scattering. Note that is a Laplace transformation in E itself, and not of ^ 
like in the diffusive limit. 
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